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ABSTRACT 

There is strong evidence that the period-luminosity (PL) relation for the Large 
Magellanic Cloud (LMC) Cepheids shows a break at a period around 10 days. Since the 
LMC PL relation is extensively used in distance scale studies, the non-linearity of the 
LMC PL relation may affect the results based on this LMC calibrated relation. In this 
paper we show that this problem can be remedied by using the Wesenheit function in 
obtaining Cepheid distances. This is because the Wesenheit function is linear although 
recent data suggests that the PL and the period-colour (PC) relations that make up 
the Wesenheit function are not. We test the linearity of the Wesenheit function and 
find strong evidence that the LMC Wesenheit function is indeed linear. This is because 
the non-linearity of the PL and PC relations cancel out when the Wesenheit function 
is constructed. We discuss this result in the context of distance scale applications. We 
also compare the distance moduli obtained from //o = Hv — R(^v — Hi) (equivalent to 
Wesenheit functions) constructed with the linear and the broken LMC PL relations, 
and find that the typical difference in distance moduli is ~ ±0.03mag. Hence, the 
broken LMC PL relation does not seriously affect current distance scale applications. 
We also discuss the random error calculated with equation /xo = Hv — R(Hv — Hi), and 
show that there is a correlation term that exists from the calculation of the random 
error. The calculated random error will be larger if this correlation term is ignored. 
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1 INTRODUCTION 

The Cepheid period-luminosity (PL) relation plays a major 
role in distance scale studies, which can ultimately be used 
to determine the Hubble constant. The calibrating PL rela- 
tion currently used is based mainly on the Large Magellanic 
Cloud (LMC) Cepheids, as applied by the Ho Key Project 
team |F reedman et al.|2 00l|) as well a s in other studies (e.g., 
ISaha et alJl200lt iKanbur et al.ll2005tl . The Cepheid PL re- 
lation has long been considered to be a linear function of 
log(P) within the range of log(P) ~ 0.3 to log(P) ~ 2.0, 
where P is the pulsation period in days. 

However, the non-linea rity of the LM C PL rela- 
tion has been proposed by iTammann et alJ ll2002l) and 
IKanbur fc Ngeowl i2004l) . i.e. the LMC data are more con- 
sistent with two PL relations and a discontinuity at a period 
around 10 days. This is illustrated in the lower panels of Fig- 
urc^for the extinction corrected F-band LMC PL relation. 
The existence of two LMC PL relations is further supported 
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by the results from a rigorous statistica l test (the F-test), 
as presented in IKanbur fc Ngeowl i2004l) . which shows that 
the V- and /-band LMC PL relations are better described 
by tw o P L relations. Since the w ork of ITammann et alJ 
J2002I) and lKanbur fc Ngeowl 12004) are based on the OGLE 
(Optical Gravitationa l Lensing Experiment) LMC Cepheids 
llJdalski et alJll999bl) . which are truncated at lo g(P) ~ 1.5, 
ISandage et al. i2004) and IKanbur et alJ <l2005l) used addi- 
tional data that are available from the literature, especially 
those with log(P) > 1.5, to further support the existence 
of two PL relations in the LMC. These studies of the non- 
linear LMC PL relations are focused on the optical bands, 
as they are mainly based on the OGLE data, we hence dis- 
cuss the non-linear LMC PL relation in the optical bands in 
this paper. The non-linear V-band LMC PL relation as seen 
from the OG LE data has been verified with the MACHO 
y-ban d data ijNgeow et al.ll200. t tl . In addition. iNgeow et alJ 
( 2005) also extended the study of the non- linear LMC PL 
relations to the JHK-bands with the 2MASS data. In Ta- 
bles we collect the slopes and the zero-points (ZP) 
for the long (log[P] > 1.0) and short period optical PL rela- 
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Figure 1. Comparison of the Wesenheit function with the extinction corrected F-band PL relation for the LMC Cepheids. On the top 
panels, we plot the Wesenheit function for the LMC data, which clearly show that the Wesenheit function is linear. We fit the data 
with a single line regression, using the least-squares method, as indicated on the figures. On the bottom panels, we show the extinction 
corrected F-band PL relations, and separate the data into the short (open circles) and long period (filled circles) Cepheids. The fitted 
PL relations for long, short and all period Cepheids are drawn as dashed, solid and dotted lines, respectively. We extend the short period 
PL relatio n to the longer period ran ges i n order to compare w ith the long period PL relations. The data in the left and right panels are 
taken from lKanbur Ngeowl 120041) and lKanbur et all l200f]t) . respectively. 



tions that are available in the literature 1 , respectively. For 
completeness, we also include the linear, unbroken PL rela- 
tion obtained using all Cepheids in these tables. The slopes 
and ZPs from different papers are consistent and agree with 
each other. Note that these LMC PL relations have been 
corrected for reddening. The reason that the LMC PL re- 
lation is non-linear is because the period-colour (PC) re- 
lation for L MC Cepheids is also non-linear across the 10 
days period jT ammann et al.||2002||Ka nbur fc Ngeowll2004t 
ISandaee et al J 12004 iKanbur et alll2005l) . The detailed in- 
vestigation of the physics behind the broken LMC PL and 
PC relation is beyond the scope of this paper, but it is of 
great interest for the studies of stellar pulsation and evolu- 
tion. 

The existence of two LMC PL relations suggests that 
in future distance scale studies, the appropriate LMC PL 
relation may need to be applied to t he long and short pe - 
riod Cepheids, respectively (see, e.g.. IKanbur et aT1l2003f> . 
However, all of the previous applications of the LMC PL 
relation were based on the linear version (in a sense, it is 
an approximation of two PL relations). Hence an immedi- 
ate question that arises is: how does the existence of two 
LMC PL relations affe ct previous s tudies? (A similar ques- 
tion was also asked bv lFeastll2003l ) In this paper we show 



1 The /-band PL relations from lKanbur fc Ngeowl <2004h are fit- 
ted with slightly different definition of the /-band mean mag- 
nitudes. Hence we refit the /-band PL relations with the con- 
ventional and reddening corrected /-band mean magnitudes. The 
results from the F-test remain unchanged. 



that if the Cepheid distance to a galaxy is derived using the 
Wesenheit function (or the equiv alent lip in equation [41 ), for 
example in the Ho Key Project l|Freedman et alJl200lL and 
references therein) or in the A raucaria Project (IGieren et alJ 
12004 IPietrzvriski et al ]|2004l) . then the results might not be 
affected (see Section 2.3 for details). This is because the 
Wesenheit function is linear, even though the PL and PC 
relations that make up the Wesenheit function are not, as 
shown in Section 2.2. Given that the LMC PL relation is not 
linear, and the Wesenheit function has been applied in many 
papers, we would like to examine the effect of broken LMC 
PL and PC relations on the application of the Wesenheit 
function to the distance scale. Again, this will be examined 
in the optical bands because the V- and /-band LMC PL 
relation s have been frequently applied in the literature (as 
in, e.g.. lFreedman et alJl200Jl . 

In addition to the non-linearity of the LMC PL 
relation, there are some recent studies also suggest- 
ing that the Cepheid PL relation is not universal, i.e., 
the Galactic PL relation is steeper than the LMC 



counterparts 


( Tammann et alJ 12001 1 


7 ouaue et al.l 12003 ; 


Kanbur et alJ 


2003 




Ngeow & Kanbur 


2004 




Storm et alJ 


2004). It is possible that the Wesenheit 


Function may 



or may not depend on metallicity: Some preliminary 
studies for bo th viewpoints can be f ound in the litera- 
ture (see, e.g.. iMoffett fc Barnedll98& ICaputo et al]l2000): 
[ Baraffe fc AlibertfeOOltlPietrzvnski et all2004IStorm et alJ 
2004). However the detailed study of the metallicity depen- 
dency of the Wesenheit function is beyond the scope of this 
paper. This paper only studies the linearity of the Wesenheit 
Function for the LMC Cepheids. 



LMC Wesenheit Function 

Table 1. Slopes for the LMC PL relations 51 . Long and short periods are referred to P > 10 and P < 10 days, respectively. 
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Band 


Period-Range 


TR02 


STR04 


KN04 


KNB05 


V 


All 


-2.760 ±0.031 b 


-2.702 ±0.028 


-2.746 ±0.043 


-2.736 ±0.036 


V 


Short 


-2.86 ±0.05 


-2.963 ± 0.056 


-2.948 ± 0.065 


-2.937 ±0.060 


V 


Long 


-2.48 ±0.17 


-2.567 ±0.102 


-2.350 ±0.252 


-2.598 ±0.161 


I 


All 


-2.962 ±0.021 b 


-2.949 ± 0.020 


-2.965 ± 0.028 


-2.965 ± 0.024 


I 


Short 


-3.03 ±0.03 


-3.099 ± 0.038 


-3.096 ± 0.043 


-3.090 ± 0.041 


I 


Long 


-2.82 ±0.13 


-2.822 ± 0.084 


-2.737 ±0.179 


-2.918 ±0.112 


B 


All 




-2.340 ± 0.037 






B 


Short 


-2.42 ±0.08 


-2.683 ± 0.077 






13 


Long 


-1.89 ± 0.62 c 


-2.151 ±0.134 







a The references are: TR02 = lTammann fc ReindJ <2002l) ; STR04 = ISandage et alj |20oi); KN04 = lKanbur fc Ngeowl J2004I) : KNB05 
= lKanbur et al.U2005l) 

b Sinc e [janimann & Rcindl (2002) docs not give the results from the fit to all Cepheids, we adopted the slopes from ludalski et alj 
il999al) because the same dataset is used in both papers. 
c The number of long period Cepheids in B-band is ~ 13, hence the error is larger than the others. 

Table 2. Same as Table HI but for the zero-points of the LMC PL relations, by assuming fiLMC = 18.50mag. 



Band 


Period-Range 


TR02 


STR04 


KN04 


KNB05 


V 


All 


-1.458 ±0.021 a 


-1.451 ±0.022 


-1.401 ±0.030 


-1.412 ± 0.025 


V 


Short 


-1.40 ±0.03 


-1.295 ±0.036 


-1.284 ±0.041 


-1.295 ± 0.038 


V 


Long 


-1.75 ±0.20 


-1.594 ±0.135 


-1.795 ±0.298 


-1.523 ± 0.198 


1 


All 


-1.942 ±0.014 a 


-1.896 ±0.015 


-1.889 ±0.019 


-1.890 ±0.017 


1 


Short 


-1.90 ±0.02 


-1.806 ± 0.024 


-1.813 ±0.027 


-1.817 ±0.026 


I 


Long 


-2.09 ±0.15 


-2.044 ±0.111 


-2.109 ±0.212 


-1.909 ± 0.138 


B 


All 




-1.160 ±0.029 






B 


Short 


-1.18 ±0.05 


-0.955 ±0.049 






B 


Long 


-1.65±0.74 b 


-1.364 ±0.177 







a Since Tammann & Reindl 1 2002) does not give the results from the fit to all Cepheids, we adopted the zero-points from lUdalski et alj 
il999al) because the same dataset is used in both papers. 
b The number of long period Cepheids in B-band is ~ 13, hence the error is larger than the others. 



2 THE WESENHEIT FUNCTION AND ITS 
APPLICATION IN DISTANCE SCALE 

2.1 Definition 

The Wesenheit func tion (e.g., see iFreedmar] 



Freedniane^il]^^^^J& 

19821: IMadore fc Freedmanl 119911 IMoffett fc Barne; 



Tanvirl 1199' 



1988; 



1976, 
1986; 



lUdalski et all Il999al) is defined as W 
magnitude — R x colour, where R is the ratio of total-to- 
selective absorption that has to be adopted. Note that the 
definition of W depends on the adopted R which may be dif- 
ferent in different environments. A few variations of W used 
in the literature with different combinations of magnitudes 
and colours are: 



wY 



V - R(B -V), R — A V /E(B - V), (1) 

V - R(V - I), R = A V /E(V - /), (2) 
I~R(V -I), R = Ri, (3) 



where B, V & 7 denote the (intensity) mean magnitudes. 
Similar definiti ons of W in near infra red bands can be found, 
for example, in Perssoii^^^Jj200j). The other definition of 
W, as given bv Ivan den Berghl dl975l) . is different than the 
one given in equation (1). The van den Bergh version of W 
replaces R by the slope of the constant-period line in the 



colour - magnit u de diagram (CMP'). [ Madore fc FreedmarJ 
J.1991F) (also in IMoffett fc Barnesl Il986t) have pointed out 
some problems with the van den Bergh version and the ad- 
vantage of using R in the definition of W . The biggest ad- 
vantag e of using W is that it is r eddening-free (see, for ex- 
ample. lMadore fc FreedmarJ 1 99 ll) . i.e. W = V-R(B-V) = 
V -R(B-V)o = Wo, where V a and (B-V) denote the in- 
trinsic visual magnitude and colour, as the effect of interstel- 
lar extinction on the observed magnitude and colour cancel 
out for a star (not only for Cepheids). Another advantage of 



using W is that the scatter in the W-\or(P) 


plot is reduced 


jMadord Il982l: IMadore & FreedmarJ 


19911: 


Bohm-Vitensd 


1997llTanvirlll997l.ll999l: lUdalski et al. 


1999a 


), as compared 



to the scatter in the V- or /-band PL relations (see Figure^. 
The remaining scatter is due to the combination of photo- 
metric errors and the finite width of the instability strip (for 
example, see iBrodie fc Madorelll98d: IMadore fc FreedmarJ 
ll99ltlB5hm-Vitensell997l : lGieren et al.ll998t) . Furthermore, 
an equivalent definition of W with the combination of ab- 
solute magnitude and colour can be formulated, as Wm = 
My — R(My — Mi), then the distance modulus can be ob- 
tained, i.e. nw — W — Wm- It is straight forward to show 
that the following equation: 



/Jo 



fiv - R(p.v - A*/) 



(4) 
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Table 3. Comparison of the slopes for W obtained using the 
Cepheids with short, long and all periods. 



Dataset 


Short Long 


All 


V-(V - I) Combination, R = 2.45. 


TR02 
STR04 
KN04 
KNB05 


-3.276 ±0.103 -3.313 ±0.403 - 
-3.296 ±0.124 -3.192 ±0.253 - 
-3.311 ±0.141 -3.298 ±0.571 - 
-3.312 ± 0.133 -3.382 ±0.360 - 


3.255 ± 0.068 
3.307 ±0.064 
3.283 ± 0.093 
3.297 ± 0.079 


I-(V - I) Combination, R = 1.55. 


TR02 
STR04 
KN04 
KNB05 


-3.293 ±0.109 -3.347 ±0.423 - 
-3.310 ±0.130 -3.217 ± 0.266 - 
-3.325 ±0.149 -3.337 ±0.601 - 
-3.327 ±0.140 -3.414 ±0.379 - 


3.275 ± 0.072 
3.332 ± 0.067 
3.304 ± 0.098 
3.320 ±0.083 


V-(B - V) Combination, R = 3.24. 


TR02 
STR04 


-4.286 ±0.335 -4.392 ±2.134 
-3.870 ±0.344 -3.915 ± 0.613 - 


3.875 ± 0.169 



Figure 2. The residuals of W from Figure HI The dashed line 
indicates the zero residuals. The residuals are scattered around 
the dashed lines, indicating that the Wesenheit function is linear. 

is equivalent to [iw, where equation (4) is frequently applied 
in determini ng the extra-galactic Cepheid distances (see, 
for e xample. lAllen fc Shanks! 120041; iFreedman et al.l l200li 



iKanbur et alJ 120031 ISaha et all l200ll ; iTanvir et al.lll999l) 

Therefore, using equation (4) to obtain the distance modu- 
lus is equivalent to obtaining the distance modulus by fitting 
the W-log(P) plane with the empirical VFM-log(P) relation. 



2.2 Testing the Linearity of the Wesenheit 
Function 

Both the Wesenheit function and equation (4) can be writ- 
ten as a combination of V- and /-band PL relations, and 
they are adopted from the LMC PL relations to derive dis- 
tances. However, as we mentioned in the Introduction, the 
LMC PL relations are not linear, hence the applicability of 
the Wesenheit function and equation (4) is immediately in 
question. In this sub-section we would like to test the lin- 
earity of the Wesenheit function as follows. 

The PL relation in bandpass A can be written as: M\ = 
a* + [3* log(P). The superscript X denotes the adopted 
period range, which is either for short (S, log[P] < 1.0), long 
(L, log[P] > 1.0) or all (A, short±long) periods. Then the 
V-(V — I) Wesenheit function becomes (similar expressions 
can be derived for other magnitude-colour combinations): 



W x = (1 - R)av + Raf + [(1 - 
The linearity of W demands that 



R)8$ + R3?]log(P), (5) 



(1 
(1- 



R)$ + R/3? 



R)al 



Ra S ! 



(1 
(1 



■R)pv 

R)a^ 



Rf3f 
- Raf 



(for 
(for 



slope) 
ZP). 



(6) 
(7) 



By using the slopes in Table Q we can calculate the val- 
ues of the left-hand side and right-hand side in equation 
(6), as well as the slope for W A by using the unbroken 



+ R^ajj . The 



PL relation with equation (5). The results are summa- 
rized in Table |3 with different magnitude-colour combina- 
tions. The adopted values for R in these combinatio ns are: 
R = A V /E(B - V) = 3.24 JUdalski et al.lll999bl): R = 
Ay/E (V-I) = 2.45 iFreedman et al.l200ll:lTanvir fc Boyle! 
Il999lk and R = Ri = 1.55 JUdalski et al.lll999al) . The er- 
rors in Tableware estimated with the standard formula for 
propagation of errors, i.e. a 2 — (1 - 
same is done for the ZP in Table ^] 

It can be seen immediately from Table |H] & 0] that the 
short period Wesenheit function is consistent with the long 
period Wesenheit function, as demanded by equation (6) 
& (7). Therefore, the Wesenheit function can be regarded 
as a linear function of log(P). Furthermore, the short and 
long period Wesenheit functions are also consistent with the 
Wesenheit function obtained from using all Cepheids in the 
LMC or the linear, unbroken PL relation. Note that the 
value of ~ —3.3 for the slope of the Wesenheit function (fiw) 
w ith I-(V-I) combinat ion also agrees with the val ues given 
inlUdalski et all (Il999al. f3 w = -3.277±0.014) or in lUdalskil 
fcOOd. Bw = -3.300±0.011'). The linearity of the Wesenheit 
function can be immediately seen from the top panels of 
Figure The residuals of the W from the fitted regressions 
in Figure are also plotted out as function of period in 
Figure |5] If the Wesenheit function is non-linear and can 
be broken into long and short period Wesenheit functions, 
as in the LMC PL or PC relations, then the residual plots 
are expected to show a trend for the long perio d Cepheids 
(as in the figure 4 from IKanbur fc Ngeowlf2004l for the PC 
relation). However, there is no obvious trend of the residuals 
seen from Figurc|2] This further supports the linearity of the 
Wesenheit function. 

A better and more sophisticated test of the linearity i s 
using the P-test iWeisberdll98ol : IKanbur fc Ngeowl I2004F) . 
The null hypothesis in our P-test is that a single linear re- 
gression is sufficient, while the alternate hypothesis is that 
two linear regressions are needed to describe the data. These 
regressions can be obtained with the standard least squares 
regression met hod. The setup and the f ormalism for the P- 
test is given in IKanbur fc Ngeowl j2004t) . By calculating the 
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Table 4. Same as Table but for the zero-points. 



Dataset 


Short Long 


All 


V-(V - I) Combination, R = 2.45. 


TR02 
STR04 
KN04 
KNB05 


-2.625 ± 0.066 -2.583 ± 0.468 - 
-2.547 ±0.079 -2.696 ± 0.335 - 
-2.580 ±0.089 -2.564 ±0.676 - 
-2.574 ± 0.084 -2.469 ± 0.444 - 


2.644 ± 0.046 
2.541 ± 0.049 
2.597 ±0.064 
2.583 ±0.055 


I-(V - I) Combination, R = 1.55. 


TR02 
STR04 
KN04 
KNB05 


-2.675 ±0.069 -2.617 ± 0.492 - 
-2.598 ± 0.083 -2.741 ± 0.352 - 
-2.633 ± 0.094 -2.596 ± 0.711 - 
-2.626 ±0.089 -2.507 ± 0.467 - 


2.692 ±0.048 
2.586 ±0.051 
2.645 ±0.067 
2.631 ±0.058 


V-(B - V) Combination, R = 3.30. 


TR02 
STR04 


-2.113 ±0.206 -2.074 ±2.543 

-2.397 ± 0.220 -2.339 ± 0.810 2.394 ±0.132 



broken PL relations. The non-linearities for both of the PL 
and PC relations cancel out when the Wesenheit function is 
formulated. 

In the next subsection, the result from the corollary 
A (equation [8]) are applied to the distance scale measure- 
ments, as well as investigate the accuracy of using the linear 
LMC PL relation vs. the broken LMC PL relation. 



2.3 Distance Scale Application 

For an ensemble of Cepheids in a target galaxy, with N 
of them being used in determining the distance, then the 
distance modulus for the i th Cepheid in bandpass A (usu- 
ally V and i) is fJ,\ = m\ — (3\ log(Pi) — ct\. The distance 
modulus in bandpass A for these Cepheids can be obtained 
by taking the unweighted mean to individual Cepheids, i.e. 
TJx = Ma- The reddening- free distance modulus for 

the i th Cepheid can be calculated with equation (4), and it 
is straight forward to show that: 



F value with N data points (note that F is a function of 
N), we can obtain the probability, p{F), under the null hy- 
pothesis, for the significance of the F value. In general, the 
null hypothesis can be rejected if F is large (e.g., F > 4) or 
p(F) is small. For example, p{F) < 0.05 indicates that the 
null hypothesis can b e rejected at 2a level. By using the 634 
Cepheids as given in lKanbur fc Ngeowl J20Q4), we found: 

F = 1.473, p(F) = 0.230 for W = V - 2.45(V - /), 
F = 0.835, p(F) = 0.432 for W = I - 1.55{V - /). 

Similarly, for the 636 Cepheids in lKanbur et alj i2005T) the 
results from the f -test are: 

F = 1.840, p(F) = 0.160 for W = V - 2A5(V - I), 
F = 1.753, p(F) = 0.174 for W = I - 1.55(1/ - /). 

Plots of the W = V — 2.45(1/ — I) as function of period for 
these two datasets are shown in the upper panels of Figure 
The relatively small values of F and the large values of 
p{F) indicate that the null hypothesis cannot be rejected. 
Hence the results from the F-test strongly suggest that the 
Wesenheit function is li near. The results wit h the extinction 
corrections (as given in lUdalski et afl ir999bl or the removal 
of the obvious outliers in upper panels of Figure Q and in 
Figure [5] are essentially the same. 

By combining the results obtained from this sub- 
section, we found that, statistically, the Wesenheit function 
for the LMC Cepheids is indeed linear, although the LMC 
PL and PC relations are not. Figure shows a comparison 
between the linear Wesenheit function and the broken V- 
band PL relation obtained from using two datasets. There- 
fore, for the LMC Cepheids, we conclude that: 

Corollary A: The Wesenheit function is statistically lin- 
ear, and the Wesenheit functions for short (P < lOdays), 
long and all (short±long) period Cepheids are approxi- 
mately the same, i.e., 

W s ~ W L ~ W A . (8) 

The fundamental reason that the Wesenheit function is lin- 
ear is because the PC relation, e.g. (V — I) = a + blog(P), 
is also broken for the LMC Cepheids, in addition to the 



1 N 

mo = ^yiri (9) 

= 777/ — R (JmF — JU), (10) 

where /xn = iA/ — R( i As — iA). This proce d ure (see, e.g., 
F^e^dmMi^t^JJ^OOlJ; iLeonard et al.l 120031 : iKanbur et alJ 
20031 lTanvirlll997ft . i.e. calculating the distance modulus for 
individual Cepheids and taking the unweighted average to be 
the final distance modulus to that target galaxy, is equiva- 
lent to fitting the V- and /-band PL relations to the data and 
obtaining the /iy and that apply to equation (4). If the 
PL relation is linear, then the single, unbroken PL relation 
(in both V- and /-band) can be adopted to fit to all Cepheids 
to obtain the distance modulus (/iff ). However, if the PL re- 
lation is non-linear, as in the case of LMC Cepheids, then the 
long period PL relation should be applied to the long period 
Cepheids, and similarly for the short period Cepheids. For 
example, attempts to use the broken LMC PL relations to 
ca librate the ext r a-gala c tic Cepheid dis t ances c an be found 
in IKanbur et alJ (1200311 : ILeonard et"aH feOOSt : iThim et alJ 
J2003ll2004h ~ 

As mentioned before, /to from equation (4) is equivalent 
to fj-w = W — Wm, i.e., /to = Hw- Since the Wesenheit func- 
tion (both W and Wm) for LMC Cepheids is linear (from 
corollary A or equation [8]), then the /to obtained from us- 
ing either the broken LMC PL relation (for long and short 
period Cepheids respectively) or the linear LMC PL rela- 
tion (for all Cepheids) is expected to agree well with each 
other. In other words, the difference in the distance mod- 
ulus (A/to) when using either the linear or the broken PL 
relation should be small. For example, the difference in dis- 
tance modulus (A/io = Mo 4 "Mo) when using either the linear 
or the broken long period LMC PL relation can be quanti- 
tatively estimated if the mean log-period (log[P]t) of the 
long period Cepheids in the target galaxy is known. For the 
distance mod ulus obtained with equation (4) & (9), A/to is 
expressed as (IKanbur et al.ll2003Tl : 

A/to = -A/3 w log(P) L - Aa w , (11) 
where Af3w and Aaw are the differences in slopes and ZPs 
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for the linear and the broken long period Wesenheit function. 
These coefficients can be calculated using the slopes and ZPs 
for the V-(V — I) Wesenheit function as given in Table|2]& 
[I] respectively. If log(P)_t ~ 1.4, then the value of A^to esti- 
mated from equation (11) is around ±0.03mag., correspond- 
ing to ~ 1.5% change in distance (in Mpc). Furthermore, for 
1.0 < log(P)i, < 2.0, equation (11) suggests that the maxi- 
mum A no when using the linear or the broken long period 
PL relation from Table is ~ ±0.07mag. A similar es- 
timation can be done for the short period Cepheids (with 
0.0 < log [P]s < 1.0) with maximum Afj, of ~ ±0.02mag. 
Even though the value of A/Xq is small, it has a systematic 
effect on the distance scale because of the different PL rela- 
tions used. 

We give two examples to illustrate the above discussion 
further: 

(i) iKanbur et alJ <l2003t) calculated the distance moduli 
to 25 HST observed galaxies with the linear and the broken 
long period LMC PL relations (as ~ 98% of the Cepheids 
detected in these galaxies have period longer that 10 days) . 
By comparing the distance moduli derived from the lin- 
ear PL relation and the broken long period PL relation, 
the unweighted average of Afio in these 25 galaxies is 
—0.021 ± 0.002mag., in the sense that the linear PL rela- 
tion systematically produces slightly smaller distance mod- 
uli than the broken long period PL relation. S ince the mean 
log-p eriod for most of these galaxies is ~ 1.4 iKanbur et alJ 
2003), then the difference of — 0.02mag. is expected and 
agrees well with the above discussion. 

(ii) To compare the distance moduli obtained from using 
the linear and the broken PL rela tions, we use the Ce pheids 
in IC 4182 (with the data from ICibson et al-lEoOol as an 
example, because this galaxy contains roughly equal number 
of short (N — 13) and long (N — 15) period Cepheids. 
The distance moduli from using the linear, unbroken PL 
relation with all 28 Cepheids (fj-^) are compared with the 
mean distance moduli (Jto, from equation [9]) obtained from 
applying the broken long and short period PL relation to 
the individual long and short period Cepheids, respectively. 
The PL relations used, include the linear and the broken 
PL relations, are taken from Table & [5] and the results 
are presented in Table |5] From this table, Ho shows a good 
agreement to JIq, with a difference of ~ O.Olmag. or smaller 2 . 

In short, when using the LMC PL relations to derive 
the Cepheid distances, we conclude that: 

Corollary B: The distance moduli are approximately 
the same when using either the linear PL relation (to all 
Cepheids) or the means with broken PL relation (to long 
and short period Cepheids respectively), i.e., 

Ho ~ Mo- (12) 

Both of the approaches will give similar and consistent dis- 
tance moduli. Recall that the linear PL relation is the ap- 
proximation of the broken PL relation, and the accuracy of 

2 Note that the median values of ~ 28.25mag. for fig and /To 
in Table 151 agree well to the distance modulus obtain ed from the 
TRGB (tip of red giant branch) method as given in lSakai et alJ 
J2004. [itrgb = 28.25 ± 0.06mag.). 



Table 5. The distance modulus (u ) to IC 4182. 



LMC PL relation 




m> 


TR02 


28.267 ±0.053 


28.277 ±0.053 


SRT04 


28.222 ±0.053 


28.217 ±0.053 


KN04 


28.250 ±0.053 


28.247 ±0.053 


KNB05 


28.253 ±0.053 


28.257 ±0.054 



using the linear PL relation is roughly ±0.03mag. Therefore, 
researchers have the freedom to use either the linear LMC 
PL relation as an approximation or using the correct broken 
LMC PL relation in deriving the Cepheid distances. 

2.3.1 A Note on the Random Errors for [iq from 
Equation (4) 

By definition, the random error for 71 = -L /it*, due to 

purely statistical fluctuations, is: 

JV 
i — l 

This equation holds for the JT\ and /ZTj (from equation [9]). 
We can expand the expression for a 2 by substituting the 
expression for /i and /lo to equation (13), then: 

o-o = (1- Rf°v + R 2 o-j -CORRvi, (14) 

where Oy and a 2 are given by equation (13), and CORRv/ = 
2 n(n-i) ~~Pv){lA ~~ W)i a term for the correlated 

resid uals from both bands jTanviJ Il997t iFreedman et alJ 
2001). Note that the use of equation (13) to estimate the 
random errors for the Cepheid dista nce modulus has been 
practiced in the lit erature (e.g., see IFreedman et al.l l200ll : 
IKanbur et afll2003lL 

On the other hand, if we apply the standard equa- 
tion for the propagation of errors (POE) to equation (4) 
or (10), by ignoring the correlation term, then we have 
ao(POE) = (1 - R) 2 cf v + R 2 aj. By comparing this expres- 
sion to equation (14), we obtain: a 2 (POE) = ct 2 + CORR V j. 
One can immediately see that a 2 {POE) is greater than er 2 
since the term CORRv/ is mostly likely to be positive. 
This is because (a) ii > 1 from the extinction curve; and 
(b) y^(/Hy — JTv)(h\ ~ ~pi) is mostly likely to be positive. 
The second condition is due to the existence of the period- 
luminosity-colour relation: if the V-band magnitude for a 
Cepheid is above (or below) the ridge-line of the F-band 
PL relation, then the corresponding /-band magnitude will 
also be above (or below) the ridge-line of the /-band PL 
relation. The result is that if (fj,y ~ T^v) ls positive (or nega- 
tive), then —~pi) is also positive (or negative), and hence 
the product of these two is positive. Further, if the extinc- 
tion and/or the correlated errors (from measurement) make 
(fly — ~pv) to be positive/negative, then (p\ — JTJ) is also 
going to be positive/negative, and again the product will 
be positive. Errors estimated from error propagation with 
equation (10), i.e. ao(POE), will ignore the CORRvi term 
and hence resulted a larger random error than the random 
error estimated from equation (13). Note that all the errors 
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discussed here, the a, are random errors only, which do not 
include the systematic errors such as the errors arise from 
the calibrations, the width of the instability s trip, and others 
(see, e.g., the discussion bv lSaha et alJ I^OOCO. 



3 CONCLUSION 

Due to the recent disc overy of the n on-linearity for the LMC 
PL and PC relations, iFeaslJ l)2003l) has asked the following 
critical question: 

"■ ■ • is there a significant slope difference between short and long 
(>~ lOdays) Cephcids that would seriously affect the calibration 
and use of PL [and PC] relation?" 

From this study, we showed that this problem can be reme- 
died with the application of the Wesenheit function in dis- 
tance scale studies. This is because the Wesenheit function 
for the LMC Cepheids is linear (corollary A), as shown in 
Section 2.2, although the LMC PL and PC relations are 
not. Therefore, the Cepheid distances obtained with the 
Wesenheit function or the equivalent fio would not be af- 
fected with the recent finding of the broken PL and PC 
relations. We also found that the typical difference in dis- 
tance modulus from using the linear or the broken PL rela- 
tions is about ±0.03mag. Hence, researchers can choose to 
apply either the linear or the broken LMC PL relations to 
obtain the Cepheid distances, without worrying that these 
two approaches will give inconsistent results, as both ap- 
proaches are equally applicable in deriving the Cepheid dis- 
tance (corollary B). Therefore, the broken PL relation found 
in the LMC Cepheids will not seriously affect the previous 
applications of the linear LMC PL relation in distance scale 
st udies because the effect is minimal. The question raised 
bv lFeastl l|2003t) is essentially answered. 



ACKNOWLEDGEMENTS 

We would like to thank the useful discussion with W. Gieren, 
D. Leonard, S. Nikolaev and G. Tammann. We would also 
like to thank the anonymous referee for useful suggestions. 



REFERENCES 

Allen, P. & Shanks, T., 2004, MNRAS, 347, 1011 
Baraffe, I. & Alibert, Y., 2001, A&A, 371, 592 
Bohm-Vitense, E., 1997, AJ, 113, 13 
Brodie, J. & Madore, B., 1980, MNRAS, 191, 841 
Caputo, F., Marconi, M. & Musella, I., 2000, A&A, 354, 
610 

Feast, M., 2003, in Stellar Candles for Extra- Galactic Dis- 
tance Scale, ed. D. Alloin & W. Gieren, Lecture Notes in 
Physics Vol. 635, Springer, New York, pg. 45 

Fouque, P., Storm, J. & Gieren, W., 2003, in Stellar Can- 
dles for Extra- Galactic Distance Scale, ed. D. Alloin & W. 
Gieren, Lecture Notes in Physics Vol. 635, Springer, New 
York, pg. 21 

Freedman, W., 1988, ApJ, 326, 691 

Freedman, W., Wilson, C. & Madore, B., 1991, ApJ, 372, 
455 



Freedman, W., Madore, B., Gibson, B., Ferrarese, L, Kel- 
son, D., Sakai, S., Mould, J., Kennicutt, R., et al., 2001, 
ApJ, 553, 47 

Gibson, B., Stetson, P., Freedman, W., Mould, J., Kenni- 
cutt, R., Huchra, J., Sakai, S., Graham, J., et al., 2000, 
ApJ, 529, 723 

Gieren, W., Fouque, P. & Gomez, M., 1998, ApJ, 496, 17 
Gieren, W., Pietrzyhski, G, Walker, A., Bresolin, F., Min- 

niti, D., Kudritzki, R.-P., Udalski, A., Soszyhski, I., et al., 

2004, AJ, 128, 1167 
Groenewegen, M. A. T., 2000, A&A, 363, 901 
Kanbur, S., Ngeow, C, Nikolaev, S., Tanvir, N. & Hendry, 

M., 2003, A&A, 411, 361 
Kanbur, S. & Ngeow, C, 2004, MNRAS, 350, 962 
Kanbur, S., Ngeow, C. & Buchler, R., 2005, in-preparation 
Leonard, D., Kanbur, S., Ngeow, C. & Tanvir, N., 2003, 

ApJ, 594, 247 

Madore, B., 1976, Royal Greenwich Obs. Bull. 182, 153 

Madore, B., 1982, ApJ, 253, 575 

Madore, B. & Freedman, W., 1991, PASP, 103, 933 

Moffett, T. & Barnes, T., 1986, ApJ, 304, 607 

Ngeow, C. & Kanbur, S., 2004, MNRAS, 349, 1130 

Ngeow, C, Kanbur, S., Nikolaev, S., Buonaccorsi, J., Cook, 

K. & Welch, D., 2005, MNRAS submitted 
Persson, S., Madore, B., Krzemihski, W., Freedman., W., 

Roth, M. & Murphy, D., 2004, AJ, 128, 2239 
Pietrzyhski, G, Gieren, W., Udalski, A., Bresolin, F., Ku- 
dritzki, R.-P., Soszyhski, I., Szymahski, M. & Kubiak, M., 
2004, AJ, 128, 2815 
Saha, A., Labhardt, L. & Prosser, C, 2000, PASP, 112, 163 
Saha, A., Sandage, A., Tammann, G, Dolphin, A., Chris- 
tensen, J., Panagia, N. & Macchetto, F., 2001, ApJ, 562, 
314 

Sakai, S., Ferrarese, L., Kennicutt, R. & Saha, A., 2004, 
ApJ, 608, 42 

Sandage, A., Tammann, G. A. & Reindl, B., 2004, A&A, 
424, 43 

Storm, J., Carney, B., Gieren, W., Fouque, P., Latham, D., 

Fry, A., 2004, A&A, 415, 531 
Tammann, G. A. & Reindl, B., 2002, Ap&SS, 280, 165 
Tammann, G. A., Reindl, B., Thim, F., Saha, A. & 

Sandage, A., 2002, in A New Era in Cosmology, eds. N. 

Metcalfe & T. Shanks, ASP Conf. Series Vol. 283, Astron. 

Soc. Pac, San Francisco, pg. 258 
Tammann, G. A., Sandage, A. & Reindl, B., 2003, A&A, 

404, 423 

Tanvir, N., 1997, in The Extra-galactic Distance Scale, Eds. 
Livio, M, Donahue, M. & Panagia, N., Cambridge Univer- 
sity Press, Cambridge, pg. 91 

Tanvir, N., 1999, in Harmonizing Cosmic Distance Scales 
in a Post-HIPPARCOS Era, eds. Egret, D. & Heck, A., 
ASP Conf. Series Vol. 167, Astron. Soc. Pac, San Fran- 
cisco, pg. 84 

Tanvir, N. & Boyle, A., 1999, MNRAS, 304, 957 

Tanvir, N., Ferguson, H. & Shanks, T., 1999, MNRAS, 310, 
175 

Thim, F., Tammann, G, Saha., A., Dolphin, A., Sandage, 
A., Tolstoy, E. & Labhardt, L., 2003, ApJ, 590, 256 

Thim, F., Hoessel, J., Saha, A., Claver, J., Dolphin, A. & 
Tammann, G, 2004, AJ, 127, 2322 

Udalski, A., Szymahski, M., Kubiak, M., Pietrzyhski, G, 
Soszyhski, I., Wozniak, P. & Zebrun, K, 1999a, AcA, 49, 



8 Ngeow & Kanbur 



201 

Udalski, A., Soszyriski, I., Szymariski, M., Kubiak, M., 
Pietrzyhski, G., Wozniak, P. & Zebrun, K., 1999b, AcA, 
49, 223 

Udalski, A., 2000, AcA, 50, 279 

van den Bergh, S., 1975, Stars & Stellar System Vol. IX, 
Galaxies and the Universe, Eds. Sandage, A., Sandage, M. 
& Kristian, J., University of Chicago Press, Chicago, pg. 
509 

Weisberg, S., 1980, Applied Linear Regression, John Wiley 
& Sons, I s * Ed. 



